microxanox_version_here <- "0.9.1"
microxanox_version_path <- paste0("microxanox_v", microxanox_version_here)
if (packageVersion("microxanox") < package_version(microxanox_version_here)) {
stop(paste0("microxanox version needs to be at least ", microxanox_version_here))
}
print(paste0("Version of microxanox: ", microxanox_version_here))
## [1] "Version of microxanox: 0.9.1"
wait_time <- "1e+06"
ss_9s <- readRDS(here::here("data",
microxanox_version_path,
"temporal_method",
"event_definition_2",
paste0("ss_data_-1strains_waittime", wait_time, "_event_definition_2.RDS")))
max_CB_var_gmax <- max(ss_9s$CB_var_gmax_s)
max_SB_var_gmax <- max(ss_9s$SB_var_gmax_s)
max_PB_var_gmax <- max(ss_9s$PB_var_gmax_s)
ss_9s <- ss_9s %>%
mutate(standvar_CB_gmax = CB_var_gmax_s / max_CB_var_gmax,
standvar_SB_gmax = SB_var_gmax_s / max_SB_var_gmax,
standvar_PB_gmax = PB_var_gmax_s / max_PB_var_gmax)
this_dat <- ss_9s %>%
filter(
SB_var_gmax_s == 0,
PB_var_gmax_s == 0
)
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
CB_var_gmax_s == 0,
PB_var_gmax_s == 0
)
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
CB_var_gmax_s == 0,
SB_var_gmax_s == 0
)
this_dat <- this_dat[order(this_dat$PB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
PB_var_gmax_s == 0,
abs(standvar_CB_gmax - standvar_SB_gmax) < 0.00001
)
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
SB_var_gmax_s == 0,
abs(standvar_CB_gmax - standvar_PB_gmax) < 0.00001
)
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
CB_var_gmax_s == 0,
abs(standvar_SB_gmax - standvar_PB_gmax) < 0.00001
)
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])
this_dat <- ss_9s %>%
filter(
abs(standvar_CB_gmax - standvar_SB_gmax) < 0.00001,
abs(standvar_SB_gmax - standvar_PB_gmax) < 0.00001
)
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]
for(i in 1:20)
print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))
#fig_state_vs_o2diff(only_SB_variable[13,])